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ABSTRACT 

We compared the quantified morphological properties of 37 intermediate redshift (0.3 < z < 1) Mgll ab- 
sorption selected galaxies to the properties of the absorbing halo gas [0.03 < W r (2196) < 2.90 A], observed 
in the spectra of background quasars. The galaxy morphologies were measured using GIM2D modeling of 
Hubble Space Telescope WFPC-2 images and the absorbing gas properties were obtained from HIRES/Keck 
and UVES/VLT quasar spectra. We found a 3.1 a correlation between galaxy morphological asymmetries 
normalized by the quasar-galaxy projected separations, A/D, and the Mgll rest-frame equivalent widths. Sat- 
uration effects cause increased scatter in the relationship with increasing W r (2196). We defined a subsample 
for which the fraction of saturated pixels in the absorption profiles is f sat < 0.5. This criterion resulted in 
a subsample of 28 systems with W r (2196) < 1.39 A. The correlation strengthened to 3.3 a. We also find a 
paucity of small morphological asymmetries for galaxies selected by Mg II absorption as compared to those of 
the general population of field galaxies, as measured in the Medium Deep Survey. The K-S probability that the 
two samples are drawn from the same galaxy population is ruled out at a 99.8% confidence level. Based upon 
four different measures of galaxy asymmetry, it is evident that the morphological perturbations of galaxies 
selected by Mgll absorption are "minor" and centrally concentrated. The A/D- W r (2796) correlation suggests 
a connection between the processes that perturb galaxies and the quantity of gas in their halos, normalized by 
the impact parameter. Since the perturbations are minor, it is clear that dramatic processes or events are not 
required for a galaxy to have an extended halo; the galaxies appear "normal". We suggest that common, more 
mild processes that populate halos with gas, such as satellite galaxy merging, accretion of the local cosmic 
web, and longer-range galaxy-galaxy interactions, consequently also induce the observed minor perturbations 
in the galaxies. 

Subject headings: quasars: absorption lines — galaxies: halos — galaxies: interactions 
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1. INTRODUCTION 

Metal-line absorbing galaxies are selected by their gas 
cross sections, i.e., the presence of extended, associated 
gaseous halos detected in quasar absorption lines. Since the 
kinematic, chemical, and ionization conditions of the gaseous 
halos can be studied in detail and compared directly to the 
galaxies themselves, these galaxy-absorber pairs provide in- 
sight into the role of gas in galaxy evolution. A connection 
between galaxy morphological characteristics and absorption 
properties would provide constraints on competing scenarios 
of the dynamics, enrichment, and geometry of halos. 

The use of high spatial resolution imaging, such as WFPC- 
2/HST and ACS/HST, provides the opportunity to study the 
morphology of intermediate redshift galaxies with unprece- 
dented detail. One of the first techniques for classifying and 
quantifying the morphologies of intermediate-to- high red- 
shift galaxies is the C-A method (I Abraham et al.ll994^ . where 
C is the concentration of light and A is the asymmetry of 
the galaxies. This technique has been applied to magnitude- 
limited samples in the Hubble Dee p Field ( Abraham et al. I 
Il996ab . the Medium Deep Survey ([Abraham et al.l 11996b). 
and the Ohio State University Bright Spiral Galaxy Survey 
(IWhvte et al.ll2002h. 

ISimard et all (120021) developed GIM2D, a two-dimensional 
decomposition fitting program that models structural param- 
eters of galaxies. This program has been used to quan- 
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tify the morphologies of H ubble Deep Field North galax- 
ies dMarleau & Simard 1998), the lu minosity-size relatio n of 
field disk galaxies for 0.1 < z < 1.1 ( Simard et al.lll999b . the 
numbe r density and lu minosity function of E/S0 galaxies for 
z < 1 dim et al.ll2002|) . and the morph ologies and bar struc- 
tures of field galaxies at 0.4 < z < 1 . 1 (IZheng et al.ll2QQ5h . 

Despite the availability of these tools, there has been little 
effort put forth toward quantifying the morphology of galax- 
ies selected by metal-line absorption in quasar spectra. To our 
knowledge, there is only one study in which the morpholo- 
gies of two Qvi absorbing galaxies were quantified using 
the C-A method (Sava geTTripp. & Lulll998|) . These galaxies 
have impact parameters 100 Kj\ and 350 Kj\ kpc, respectively 
(go = 0.05, ft \ = 0). Based upon the concentration parame- 
ter, they find that these two galaxies are Sbc-Sc spirals. They 
also find that their asymmetries are lower than that of strongly 
interacting or irregular galaxies. 

Though various classes of metal-line absorbers are avail- 
able for sample selection, the Mgll AA2796,2803 dou- 
blet is arguably the best for studying the association of 
absorption with galaxies. Mgll is observable from the 
ground for intermediate redshifts, a regime where galax- 
ies can be well resolved in HST images. Mgll absorp- 
tion probes low ionization halo and is detectable in en- 
viron ments spanning five decades of Hi column densi- 
ties (Churchill et al. 2000a), from sub-Lyman limit (e.g., 
Churchill et al. 1999b) to damped Lya systems (e.g., 
iRao & Turnsh ek 2000). The association of normal, bright, 
field galaxies with Mgll absorption at m oderate impact pa- 
rameters has been well established (e.g., Bergeron & Boisse 
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ll991l:ISteidel. Dickinson. & Perssonlll994h . 

Much of our current understanding of Mgll absorp- 
tion selected galaxies is based upon the galaxy B 
and K luminosities, rest-frame B - K colors, an d im- 
pact parameters (ISteidel. Dickinson. & PerssonI Il994t 
IChurchill. Steidel. & Vod Il996h IChurchill et all l2000bh . 
The spectroscopic types of these galaxies have been 
based upon their c o lors fr o m ground-based images (e.g., 
Bergeron & Boisse |1991l: ISteidel. Dickinson. & PerssonI 
1 994; IChen & Lanzetta 200 3]), or, in a few cases, based upon 
their spec tra (e.g.,|S pinrad et al. 1993; Guillemin & Bergeron! 
119971) . ISteidel. Dickinson. & PerssonI (11994 determined 
that the average B-K color for Mgll absorbing galaxies 
is consistent with that of an Sb galaxy, though all galaxy 
colors redder than an Sd color are represented. Using Sloan 
r, and / band images toward 700 quasars, IZibetti et all 
( 2005| |2006|) also find galaxy colors similar to those of local 
Sb-c spirals based upon statistical methods. Space-based 
WFPC-2///ST images of intermediate redshift absorption 
selected galaxies qualitatively reveal that their morpho- 
logical types appear to be simil a r to t h ose of loca l spiral 
and elliptical galaxies (ISteidell fl998b IChen et al.l [200 1 ; 
Chen & Lanzetta 2003; Kacprzak, Churchill, & Steidel 2005; 
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In the local universe, where case-by-case studies can be 
very detailed, the morphologies of galaxies have provided 
important insights for understanding the characteristics of 
gaseous halos. There are several case studies where faint halo 
gas is observed in emission to a projected distance of 13 kpc; 
the majority of these galaxies have either been harassed by 
and/or are interacting with other galaxies, as is a pparent from 
their morphological asymmetries and H i warps (Puch e_etaD 
1992; S waters, Sancisi & van der Hulst 1997; Rand 2000). If 



local interacting or perturbed galaxies have more gas ejected 
into their halos, then intermediate redshift metal-line absorp- 
tion selected galaxies might likewise display some evidence 
of morphological perturbations. If so, what can a galaxy's 
morphology and/or morphological perturbations tell us about 
the absorption characteristics of the halo gas? Moreover, what 
can galaxy perturbations tell us about the mechanisms that 
populate galaxy halos with gas? 

In order to shed light on these questions, we have applied 
GIM2D on WFPC-2/ffiT images to quantify the morpholo- 
gies of known Mgll absorbing galaxies at intermediate red- 
shifts. We aim to examine whether the morphologies of 
metal-line absorption selected galaxies can provide insight 
into the quantity and distribution of gas residing in their ex- 
tended halos. We have employed several asymmetry mea- 
sures, since each invokes different methods for quantifying 
deviations from an axisymmetric, smooth morphology. Our 
study has revealed a correlation between the rest-frame equiv- 
alent width of the Mgll A2796 transition, W r (2796), and the 
galaxy morphological asymmetries normalized by the impact 
parameter, i.e., the projected galactocentric distance to the ab- 
sorption. 

This paper is organized as follows: In § [2 we present our 
sample selection methodology, the resulting sample, and the 
data reduction methods we employed. In § [3j the basic de- 
tails of quantifying the galaxy morphological asymmetries are 
presented. We present our results and analyze the differences 
between various methods of measuring morphological asym- 
metries in § HI In § [3 we discuss the implications of our 
analysis. Concluding remarks are provided in § [6] Through- 
out this paper, we adopt a Hq = 70 km s -1 Mpc -1 , Vt m = 0.3, 



and £7 a = 0.7 cosmology. 

2. GALAXY SELECTION AND DATA ANALYSIS 

Our sample of galaxies is selected by the presence of 
Mgll absorption in ground-based, optical quasar spectra, 
primarily high resolution (R = 45 , 000, FWHM ~ 6 kms" 1 ) 
HIRE S/Keck dVogt et alJl!994h or UVES/VLT dDekker et al] 
2000) data. These high signal-to-noise ratio spectra provide 
a Mgll A2796 rest-frame equivalent width detection thresh- 
old of W r (2196) ~ 0.02 A. Thus, we have been able to se- 
lect galaxies associated with a wide range of W r (2196), in- 
cluding the so called "weak" systems, i.e., W r (2796) < 0.3 A 
(IChurchill et alJH999bLl200^ . 

In the cases where we have a HIRES or UVES quasar 
spectrum of an Mgll absorber, we have measured the ab- 
sorption properties and equivalent widths from the data. The 
HIRES spectra were reduced using IRAF 4 . The UVES spec- 
tra were reduced using the standard ESO pipeline and a cus- 
tom code called the UVES Pos t-Pipeline Echelle Reduction 
(UVES POPLER. lMurphvll2006h . Analysis of the Mgll ab- 
sorption profiles was performed usi ng graphic-based inter- 
active software of our own de sign ( Churchill et al. 1999b, 
2000at IChurchill & Vogtl 1200 ll: IChurchill. Vogt. & Cha rlton 
2003) for local continuum fitting, objective feature identifi- 
cation, and for measuring absorption properties. For Mgll 
absorbers where we do not have a HIRES or UVES quasar 
spectrum, we obtained the measured WX2796) from the liter- 
ature. 

We have limited our selection of galaxies to those 
with spectroscopically measured redshifts that are consis- 
tent with the redshift of the Mgll absorption in the range 
0.3 < z < 1.0. The galaxy redshifts are either obtained 
throu g h our previous work (ISteidel. Dickinson. & PerssonI 
19941: IChurchill. Steidel. & Vogtl 119961: ISteidel et aD fl997l 



2002), or from the literature. In order to examine the galaxy 



morphologies, we further limited our sample to include only 
those galaxy-absorber pairs in quasar fields that have been 
imaged with either the F702W and/or F814W filter using 
WFPC-2/HST. The WFPC-2 images were obtained for vari- 
ous research programs and are available in the HST archive. 
The F702W filter provides a broad-band image similar to a 
rest-frame Johnson 5-band filter for galaxies at z ~ 0.6. The 
F814W filter provides an image similar to the rest-frame B- 
band for galaxies at z ~ 0.85. 

All WFPC-2/HST images were reduced using the WFPC- 
2 Associations Science Products Pipeline (WASPP 5 ). The 
WFPC-2 astromet ry is calibrated to the USN02 Catalog 
(Mon et et all [1998). Offsets for combining images are per- 
formed either by cross-correlation, jitter information, and/or 
image world coordinate system information. Robust aver- 
age imag es are obtain ed using the method of artificial skep- 
ticism (IStetsonl fl989|) . which uses a continuous weighting 
scheme derived directly from the data. WASPP data qual- 
ity verifications include photometric and astrometric accu- 
racy and correctly set zero-points. Galaxy photometry was 
performed using the Source Extractor (Sextractor) package 

4 IRAF is written and supported by the IRAF programming group at 
the National Optical Astronomy Observatories (NOAO) in Tucson, Arizona. 
NO AO is operated by the Association of Universities for Research in Astron- 
omy (AURA), Inc. under cooperative agreement with the National Science 
Foundation. 

5 Developed by the Canadian Astronomy Data Centre (CADC) 
and the Space Telescope-European Coordinating Facility (ST-ECF): 
http .-//archive, stsci. edu/hst/wfpc2/pipeline. html 
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TABLE 1 
Galaxy-Absorber Sample 



QSO Field 


^■abs 


W r (2796) Mgn a 


D (kpc) 




rriHST 


Filter 


Exp, sec 


PID 


Q0002 + 051 


0.298059 


0.246 ±0.004 1 


59.3 ±0.3 


l 


19.96 ±0.01 


F702W 


4600 


Steidel/5984 


Q0002 + 051 


0.591365 


0.102 ±0.002 1 


36.2 ±0.4 


l 


21.24 ±0.01 


F702W 


4600 


Steidel/5984 


Q0002 + 051 


0.851407 


1.119±0.013 1 


25.9 ±0.5 


l 


22.34 ±0.02 


F702W 


4600 


Steidel/5984 


QO 109 + 200 


0.5346 


2.26 2 


45.1 ±0.4 


2,3 


22.27 ±0.02 


F702W 


1800 


Disney/6303 


Q0117 + 213 


0.729075 


0.244 ±0.006 1 


55.5 =b 1.1 


1 


21.06 ±0.02 


F702W 


2008 


Zuo/6115 


QO 150 -202 


0.3887 


0.58 2 


60.0 ±0.7 


2,3 


21.14±0.01 


F702W 


5100 


Steidel/6557 


Q0229 + 131 


0.417337 


0.816 ±0.022 1 


37.5 ±0.5 


2,3 


19.72 ±0.01 


F702W 


5000 


Steidel/6557 


Q0235 + 164 


0.5240 


2.34 ±0.05 3 


12.1 ±0.6 


3 


20.30 ±0.02 


F702W 


600 


Burbidge/5096 


Q0450-132 


0.493937 


0.674 ± 0.026 1 


50.1 ±0.4 


1 


21.55±0.01 


F702W 


2500 


Steidel/5984 


Q0450-132 


0.493937 


0.674 ±0.026 1 


62.7 ±0.7 


1 


21.53±0.01 


F702W 


2500 


Steidel/5984 


Q0827 + 243 


0.524968 


2.419 ±0.017 1 


37.5 ±0.9 


4 


20.64 ±0.01 


F702W 


4600 


Steidel/5984 


Q0836 + 113 


0.7874 


2.71 4 


26.9 ±0.9 


5 


22.63 ±0.02 


F702W 


5000 


Steidel/6557 


Q1019 + 309 


0.3461 


0.70 ±0.07 5 


46.0 ±0.5 


1 


20.47 ±0.01 


F702W 


5100 


Steidel/6557 


Q1038 + 064 


0.441453 


0.673 ±0.014 1 


56.0 ±0.6 


6 


20.70 ±0.01 


F702W 


4600 


Steidel/5984 


Ql 127 -145 


0.312710 


1.773 ±0.006 1 


81.0±0.3 


2,3 


19.67 ±0.01 


F814W 


4400 


Bechtold/9173 


Ql 127 -145 


0.312710 


1.773 ±0.006 1 


45.6 ±0.3 


2,3 


19.48 ±0.01 


F814W 


4400 


Bechtold/9173 


Q1148 + 387 


0.553362 


0.644 ±0.014 1 


20.3 ±0.6 


6 


20.94 ±0.01 


F702W 


4800 


Steidel/5984 


Q1209 + 107 


0.392925 


1.210±0.006 1 


37.9 ±0.4 


3 


21.74 ±0.02 


F702W 


3600 


Bergeron/5351 


Q 1222 + 228 


0.550196 


0.098 ± 0.009 1 


36.4±1.1 


6 


22.50 ±0.02 


F702W 


5000 


Steidel/5984 


Q1241 + 176 


0.550483 


0.482 ±0.013 1 


21.4±0.4 


1 


21.40 ±0.01 


F702W 


5000 


Steidel/6557 


Q 1246 -057 


0.639908 


0.461 ±0.006 1 


29.9 ±0.8 


1 


22.21 ±0.02 


F702W 


4600 


Steidel/5984 


Q1317 + 277 


0.660053 


0.348 ±0.007 1 


103.2 ±0.6 


6 


21.35 ±0.01 


F702W 


4700 


Steidel/5984 


Q1332 + 552 


0.3740 


2.90 2 


27.9 ±0.5 


7 


19.40 ±0.01 


F702W 


2800 


Steidel/6557 


Q1511 + 103 


0.4369 


0.454 ± 0.046 6 


38.2 ±0.4 


2,3 


21.22±0.01 


F702W 


5000 


Steidel/6557 


Q1622 + 235 


0.317596 


0.491 ±0.012 1 


54.6 ±0.5 


8 


19.45 ±0.01 


F702W 


24,000 


Steidel/5304 


Q1622 + 235 


0.368113 


0.253 ±0.008 1 


114.0±0.4 


8 


20.00 ±0.01 


F702W 


24,000 


Steidel/5304 


Q1622 + 235 


0.471930 


0.769 ± 0.009 1 


34.1 ±0.5 


8 


22.27 ±0.01 


F702W 


24,000 


Steidel/5304 


Q1622 + 235 


0.656103 


1.451 ±0.011 1 


99.6 ±0.5 


8 


22.52 ±0.02 


F702W 


24,000 


Steidel/5304 


Q1622 + 235 


0.702903 


0.032 ±0.003 1 


112.7±0.9 


8 


21.63 ±0.01 


F702W 


24,000 


Steidel/5304 


Q1622 + 235 


0.797079 


0.468 ±0.009 1 


71.5±1.3 


8 


22.37 ±0.02 


F702W 


24,000 


Steidel/5304 


Q1622 + 235 


0.891280 


1.548 ±0.008 1 


23.2 ±0.5 


8 


22.65 ±0.02 


F702W 


24,000 


Steidel/5304 


Q1623 + 268 


0.8885 


0.68 ±0.03 7 


48.2 ±0.5 


1 


23.64 ±0.03 


F702W 


4600 


Steidel/5984 


Q1623 + 268 


0.8875 


0.27 ±0.02 7 


69.9 ±1.3 


1 


23.80 ±0.06 


F702W 


4600 


Steidel/5984 


Q2128-123 


0.429820 


0.464 ±0.004 1 


48.8 ±0.5 


2,3 


20.42 ±0.01 


F702W 


1800 


Mecchetto/5143 


Q2206-199 


0.751923 


0.886 ±0.003 1 


44.2 ±0.7 


3 


22.35 ±0.01 


F702W 


5000 


Steidel/6557 


Q2206-199 


0.948361 


0.253 ±0.002 1 


87.2 ±0.5 


3 


21.92 ±0.01 


F702W 


5000 


Steidel/6557 


Q2206-199 


1.017040 


1.057 ±0.005 1 


104.6 ±1.4 


3 


20.99 ±0.01 


F702W 


5000 


Steidel/6557 



a Mgll Absorption: (1) T his paper, (2) iGui llemin & Bergeron 1 1997), (3) Lanzetta & Bowen] J1992I) . (4) [Steidel, Dickinson, & Persson fl994l) 
(5) Steidel & Sargent 1 1992), (6) Foltz et al. 1 1986), (7) Sargent Boksenberg & Steidel 1 1988). b Galaxy Identification: (1) Steidel, Dickinson, & Persson 
1 1994), (2) Bergeron & Boisse 1 1991), (3) Guillemin & Bergeron 1 1997), (4) Kanekar & Chengalur 1 2001), (5) Lowenthal et al. 1 1990), (6) Steidel et al. 
1 2002), and (7) Miller, Goodrich, & Stephens 1 1987), (8) Steidel et al. 1 1997). 



(Berlin & Arnouts 1996) with a detection criterion of 1.5 a 
above background. The mfiq2w and mF%\w magnitudes wer e 
measured using the WFPC-2 zero points (IWhitmorel fl 995 ) . 
based upon the Vega system. Uncertainties in magnitudes 
range from 0.01 to 0.02. 

Our final sample comprises 37 galaxies with magni- 
tudes in the range 19.4 < mhst < 2 3.8 selected by "weak " 
(IChurchill et al.l Il999b|) "ri assic" (IChurchill et alJ l2000bh . 
and "strong" dBond et al.ll200l|) Mgll absorption spanning the 
equivalent width range 0.03 < W r (2796) < 2.90 A. Note that 
for three of the absorbers, each have two galaxies at the same 
redshifts. Since it is always possible that a pair of galaxies 
can give rise to the absorption and such conditions can pro- 
vide further insights into the galaxy-absorption connection, 
we do not exclude these pairs from our sample. 

The sample is presented in Table [TJ where we list (1) the 
quasar field, (2) the Mgll absorption redshift, (3) the rest- 
frame Mgll A2796 equivalent width, W r (2796), and uncer- 
tainty (when available), (4) the source of quoted Mgll mea- 
surements, (5) the quasar-galaxy impact parameter, D, and 
uncertainty, and (6) the reference(s) for the galaxy identifica- 
tion. We also list (7) the apparent magnitude of the galaxy and 
uncertainty as measured in the WFPC-2 filter (8), (9) the ex- 
posure time, and (10) the PI name and program identification 



of the WFPC-2 image. Further details about the quasar fields, 
the absorbing galaxies, and the absorption properties w ill be 
published in an upcoming paper dKacprzak et al .1 12007). 

The impact parameters, D, presented in Table U are com- 
puted using the galaxy 1.5 a isophotal centroids determined 
by Sextractor. There is a ~ 0.05 pixel uncertainty in the posi- 
tion of the quasar based upon centroiding errors of unresolved 
sources in the images. As described below , we modeled the 
galaxies using GIM2D ( Si mard et al.ll2002l) . which provides 
a second isophotal centroid. The dominant uncertainty in D 
is derived from the pixel offset of the galaxy isophotal cen- 
ter obtained using Sextractor and the isophotal center of the 
galaxy model determined by GIM2D. This offset is typically 
about 0.25 pixels. 

WFPC-2 portraits of the galaxies are presented in the left- 
hand panels of Figure \Tg\ The portraits have an area 10 times 
larger than the 1.5 a galaxy isophotal area. The encircled 
arrow next to the portrait provides the direction to the quasar 
relative to the galaxy; also shown are the cardinal directions. 
The accompanying panels are the image of the GIM2D model 
and the residual image. These are discussed below. 

3. QUANTIFYING GALAXY AS YMMETRIES 

For each galaxy, the morphology was quantified by fitting 
a two-component (bulge+disk) parametric model to its two- 



4 



Kacprzak et al. 



HST Image 



Residual Image 



© 



'J 



Q0002+051 z=0.298059 



e • N ;? 



Q0002+051 z=0.591365 




Q0002+051 z=0.851407 



© 



r 

^ N 



QO 109+022 z=0.5346 







Q01 17+213 z=0.729075 



© • "3 



Q0150-202 z=0.3887 



© 



Q0229+131 z=0.417337 



© 



Q0235+164 z=0.524 



© 



Q0450-132 z=0.493937 



Model Image 



Residual Image 



© 



Q0450-132 z=0.493937 







> 



Q0827+243 z=0.524968 



Q0836+113 z=0.7874 



© 



> 



Q1019+309 z=0.3461 







© 



Q1038+064 z=0.441453 







■J 



Ql 127-145 z=0.312710 



© 



Ql 127-145 z=0.312710 



© 



© 



Ql 148+387 z=0.553362 



© 



•J 



Q1209+107 z=0.392925 




Fig. la. — (left) WFPC-2 images of galaxies selected by Mgn absorption. The images are 10 times larger than the 1.5 a isophotal area. — (center) The 
GIM2D models of the galaxies, which provide quantified morphological parameters. — (right) The residual images from the models, showing the underlying 
structure and morphological perturbations of the galaxies, which are quantified using various asymmetry measures. The encircled arrow provides the direction to 
the quasar (galaxy-quasar orientation) and the cardinal directions are shown. The quasar name and redshift of Mg II absorption is shown under each set of galaxy 
WFPC-2, model and residual image. 



Mgn Absorber Galaxy Morphology 



HST Image 



Model Image Residual Image 



Q1222+228 z=0.550196 







Q1241+176 z=0.550483 







Q1246-057 z=0.639908 



Q13 17+277 z=0.660053 




Q1332+552 z=0.374 




Q151 1+103 z=0.4369 







Q1622+235 z=0.3 17596 












Q1622+235 z=0.471930 




HST Image 



Model Image Residual Image 







N E 



Q1622+235 z=0.656103 







N E 



Q 1622+235 z=0.702903 



* 



N E 



Q 1622+235 z=0.797079 







Q1622+235 z=0.891280 







E N 

V 



Q1623+268 z=0.8885 







Q1623-268 z=0.8875 




N ^ 



Q2128-123 z=0.429820 







Q2206-199 z=0.751923 




Q2206-199 z=0.948361 








Q2206-199 z= 1.0 17040 



Kacprzak et al. 



6 

dimensional surfac e brightness distribution using GIM2D 
( Sim ard et al.ll2002|) . We fitted the surface brightness of the 
bulge component with a de Vaucouleurs profile and the sur- 
face brightness of the disk component with an exponential 
profile. The centroids of the bulge and disk components are 
forced to be co-spatial; they are also fitted parameters. We 
employed this standard profile on all galaxies in our sam- 
ple in order to ensure uniform parameterization. The models 
are convolved with the WFPC-2 point spread function, which 
was determin ed at the appropriate locations on the image us- 
ing Tiny Tim (iKrist & Hookll2004l) . 

Our goal for employing GIM2D to fit HST galaxy images is 
to derive structural and morphological parameters that can be 
used for quantitative comparison with quantified absorption 
properties, primarily the Mgll rest-frame equivalent width. 
We are particularly interested in quantified measurements of 
the morphological perturbations, or asymmetries in our ab- 
sorption selected galaxies. GIM2D c omputes four measure s 
of asymmetry usin g the for malism of Abraha m et al.l (11994)) . 
ISchade et all(ll995h . and ISimard et al.l (120021) . Portraits of 
the GIM2D model galaxies (center panels) and the resid- 
ual images (right-hand panels) are presented in Figure \Ta\ 
The residual image is obtained by subtraction of the GIM2D 
model image from the HST image. For the residuals, regions 
that lie outside the Sextractor segmentation image are identi- 
cal to the HST images. 

The first asymmetry measure is obtained from the C-A 
method of lAbraham et al.l (119941) . The Abraham A is com- 
puted directly from the data (model independent) from 




ij ij 



where the 1^ are the counts at pixel ij and the 7^ 80 are the 
counts at the ij location following a 180° rotation of the 
galaxy. The domain of the sums is an elliptical isophotal area 
defined by galaxy counts 2 a above background. The back- 
ground counts, Bij, are randomly selected from pixels just be- 
yond the Sextractor segmentation image. Uncertainties in the 
A are determined by computing A for different random selec- 
tions of the B^. In GIM2D, four realizations of background 
are computed. For the best value of A, we take the average of 
the four realizations and for the uncertainty in A we compute 
the standard deviation of these realizations about the average. 

The second asy mmetry measure is the quantity Ra 
(Sch ade et al.| [l995). which is also computed using Eq.Q] ex- 
cept that it employs the GIM2D galaxy model residuals, Rfj 
and R}j°, instead of the direct counts 1^ and 7^ 80 from the 
data. The domain of the sums is also model dependent. Ra is 
evaluated for each of the GIM2D galaxy model circular aper- 
tures, 1-10 rh (half-light radii) for R t j that are 1, 3 and 5 a 
above background. Thus, there are 30 values of Ra represent- 
ing 10 circular apertures in units of for three significance 
thresholds. As with the Abraham A, the background counts 
are randomly selected from the pixels just beyond the Sex- 
tractor segmentation image. 

Two addit i onal methods, denoted as A z and D z 
(ISimard et al.l l2002h . are unique to GIM2D. Both are 
computed directly from the data, but using the GIM2D model 
circular apertures. A z is evaluated for each of ten circular 
apertures, 1-10 r^, and is the sum of all having signal 2, 3, 



and 5 a greater than that of the 180° symmetrical counterpart 
normalized by the total counts. Similar to the Ra, there are 30 
values of A z representing 10 circular apertures in units of r% 
for three significance thresholds. D z is computed exclusively 
outside the GIM2D model circular aperture 1 It is the 
sum of all 1^ that are 1.5 a above background for which 
the symmetric counterpart are less than 1.5 a above 
background. The sum is normalized by the total galaxy total 
flux. D z is a sensitive measure of highly asymmetric spiral 
arms, tidal tails, and highly asymmetric perturbations, since 
such features tend to be located to one side of a galaxy (see, 
for example, Q1246-057 and Q2206- 199 at z = 1.0170 
in Figure \T^. Comparing D z to the other asymmetric 
measures provides insight for discriminating if a galaxy has 
suffered strong interactions or only slight morphological 
perturbations. 

The WFPC-2 image exposure times range from 600 s to 
24,000 s, with a typical time of 5000 s. The quantified 
morphological parameters obtained from the 600 s exposure 
(Q0235+164) should be viewed with caution. It may also be 
of concern that the longest exposure of 24,000 s (Q1622+235) 
could have different measured morphological parameters if 
the exposure time was of the order of 5000 seconds. To verify 
that the these values did not change drastically as a function 
of exposure time, we degraded the 24,000 s exposure to both 
10,000 and 5000 s and measured the morphological param- 
eters for these three exposure times. For all seven galaxies 
in the degraded 10,000 s exposure, the measured values are 
consistent within the uncertainties of the measurement. For 
the degraded image of 5000 s, the average percent difference 
is 9% but this is not systematic (some increase slightly, some 
decrease slightly). The quoted errors in the morphological 
parameters for the Q 1622+235 field reflect the errors in the 
5000 s exposure. 

3.1. Effects of PSF Subtraction 

Seven galaxy-absorber pairs lie at relatively small im- 
pact parameters such that the quasar light potentially con- 
taminates the two-dimensional light profiles of the galaxies 
(e.g., Q0002+051 z = 0.851407, Q0235+164, Q0827+243, 
Ql 148+387, Q1241+176, Q1246-057, and Q1622+235 z = 
0.891280). In some of these cases, a diffraction spike passes 
through the galaxy or passes near the galaxy. In other cases, 
the galaxy lies within the extended tails of the complicated 
point spread function (PSF) illumination pattern of the quasar. 
Whereas galaxies far from the quasar light pattern are ob- 
served to have a uniform (flat) background, those near the 
quasar light pattern lie on a tilted and perhaps slightly non- 
uniform "background". In some cases, there is a ramping of 
the background of roughly 3-5 counts over a 60 pixel radial 
cut toward the quasar in the vicinity of the galaxy. 

In both cases, the quasar light can affect the physical bound- 
aries of the Sextractor segmentation image and/or affect the 
measured morphological parameters obtained by GIM2D. We 
examined the effects of the quasar light contamination by per- 
forming two tests. First, when possible, we modeled and sub- 
tracted the PSF of the quasar from the image and compared 
the resulting segmentation images and GIM2D model param- 
eters. Second, in cases where PSF subtraction was unsatis- 
factory, we fit a two-dimension smooth surface to the slight 
ramping of the background (due to the quasar illumination) 
around a 60 x 60 pixel region centered on the galaxy. 

In two cases (Q0827+243 and Ql 148+387), the post- and 
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pre- PSF subtracted asymmetries were consistent within the 
errors obtained via the four realizations of the background. 
In four cases (Q0002+051, Q0235+164, Q1241+176, and 
Q 1246-057), the asymmetries differed by more than their re- 
spective background induced errors (which are inherently un- 
derestimates in that Poisson noise in the galaxy counts are not 
included). However, this was not a systematic effect. In two 
of those cases (Q0002+051 and Q0235+164) the A values in- 
creased and in two cases (Q124 1+176 and Q 1246-057) the A 
values decreased. 

Due to the inherent uncertainties in PSF subtraction of sat- 
urated quasars, we adopt the A values from the pre-PSF sub- 
tracted images and use the deviations between the pre- and 
post- PSF subtracted images as a measure of the uncertainty 
in the A values. The adopted uncertainties for these seven 
galaxy-absorber pairs are roughly a factor of five larger than 
their background only induced errors. 

4. ANALYSIS AND RESULTS 

In Table [2 we present the quantified asymmetries for each 
absorption selected galaxy in our sample. In columns (1) and 
(2), we list the quasar field and the absorption redshift. In col- 
umn (3), we list the Abraham A and uncertainty. For purposes 
of later discussion, we list representative values of Ra and A z , 
for apertures 2 and for 3 a and 2 a significance thresholds, 
respectively, in columns (4)-(5) and (6)-(7). In column (8), 
we list D z . 

In order to explore possible correlations between galaxy 
asymmetries and the strength of the absorbing gas, we per- 
formed non-parametric Spearman and Kendall rank correla- 
tion tests between the various asymmetry measures, A, Ra, 
A z , and D z , and the Mgll A2796 equivalent width. We also 
examined various scalings of the asymmetries, most notably 
normalization by the impact parameter, D. 

For all tests, the null-hypothesis of no correlation could not 
be ruled out to a high level of significance (i.e., greater than 
3 a), except for the Abraham A, for which a 3.1 a correla- 
tion was found between A/D and Wr(2796). Recall that the 
Abraham A is model independent and is the only asymmetry 
based upon an elliptical aperture, whereas the Ra, A z , D z are 
derived from the circular model dependent apertures at vari- 
ous 1%. Thus, A represents a full galaxy, model independent 
asymmetry with a natural isophotal aperture. 

At face value, the correlation between A/D and W r (2196) 
for the full sample (see the inset of Figure H]) suggests that, 
once the impact parameter of absorption is normalized out, 
there may be a connection between perturbations in a galaxy's 
morphology and the quantity of Mg II absorbing gas probed 
in its halo (which is a complex distribution of column den- 
sity and velocity spread). Note that for the largest equivalent 
widths, the scatter in A/D increases substantially. Based upon 
this impression, we might speculate that there is a physically 
based transition in the behavior of W r (2196) for large equiva- 
lent widths, i.e., a transition from a "clouds in halos " scenario 
to a " galactic winds" scenario [see, for example, iBond et all 
( 2001)] . Since absorption profiles become progressively more 
saturated with increasing equivalent width, it is more conser- 
vative to simply suggest that this increase in scatter could arise 
due to a loss of information on the quantity of gas probed by 
the quasar. 

4.1. Equivalent Widths and Saturation 

The equivalent width is a single valued measure of a com- 
plex velocity distribution of gas column density. For "weak" 



and "classic" systems, distinct individual clouds are usually 
detected with low to moderate column densities, and the 
number of clouds strongly correlates with equivalent widt h 
(IChurchill & Vogtll^OOlTIChurchill. Vogt. & Charltoril2003h . 
However, for much "stronger" systems, the correlation breaks 
down for saturated profiles, indicating that the equivalent 
width can no longer be employed to accurately decipher the 
combination of column density, number of clouds, and their 
velocity distribution. 

Using a n appro a ch sim ilar to the saturation D-index em- 
ployed by lEllisonl (|2006|) to efficiently pre-selected damped 
Lya absorbers (DLAs), we computed the fraction saturated 
pixels, f sau across each HIRES or UVES Mgll A2796 ab- 
sorption profile. The full pix el range of a profil e is us ed, as 
defined using the methods of IChurchill & Vogtl (|2001). We 
use the criterion that a pixel is saturated when the counts are 
less than 1.5 (a), where (a) is the average of the uncertainty 
spectrum computed over a window ±600 km s -1 from the op- 
tical depth mean (which defines the velocity zero point). In 
Figured we present W r (2796) versus the fraction of saturated 
pixels, f sat . The solid points represent the galaxies in our sam- 
ple for which we have HIRES and UVES spectra of the Mg II 
A2796 absorption profiles. For these data, there is some sug- 
gestion that the diagram visually breaks into four quadrants 
divided by W r (2196) ^ 1.3 A and f sat ~ 0.5. These quadrants 
tend to be either fully populated or entirely void of data points. 

In order to ascertain whether this trend is universal to 
HIRES and UVES Mgll absorption profiles, we populated the 
figure with an additional 97 Mg II systems from our database 
of HIRES and UVES spectra. These systems were selected 
to cover the same redshift and equivalent width ranges of our 
samp le and will be pres ented elsewhere as part of a larger sur- 
vey (lEvans et al.l 12007). We plot these supplementary data 
as open circles in Figure [2 As can be seen, the general 
trend remains, especially the void of data in the quadrant de- 
fined by W r (2196) < 1.4 A and f sat > 0.5. The behavior of 
Mgll absorption is such that there are no profiles with a sat- 
uration fraction of f sat > 0.5 and equivalent width less than 
W r (2796) ~ 1.4 A. For equivalent widths above this value, 
there is a small population of profiles with f sat < 0.5, but none 
of these are present in our sample. 

We partitioned the ^(2796)-/^ plane into four quadrants 
in order to objectively classify the different types of absorp- 
tion profiles with the goal of defining a subsample unaffected 
by strong saturation. We refer to these as Q1-Q4, starting at 
the origin and increasing counterclockwise on the plane (Q2 is 
the empty quadrant). The final quadrants were determined by 
maximizing the number of data points in Ql and Q3 and min- 
imizing the number in Q2 and Q4. This method resulted in 
a strict equivalent width transition of W r (2196) = 1.39 A and 
a less stringent transition of saturation in the range 0.46 < 
fsat < 0.51. The dashed lines on Figure [2 define the result- 
ing quadrants, where we have adopted f sat > 0.5 to designate 
highly saturated profiles. A Kolomorov-Smirnov (K-S) test 
yielded a probability of P KS < 10~ 8 (5.2 a) that systems with 
W r (2796) < 1.39 A and W r (2196) > 1.39 A are drawn from 
the same distribution of f sat . 

There are a total of nine galaxy-absorber pairs in our sam- 
ple that fall into our "highly saturated" classification. For four 
(associated with five galaxies), we have HIRES and UVES 
spectra and we have directly computed their f sat . However, 
there are four in the range 2.3 < W r (2796) < 2.90 A for which 
we do not have HIRES or UVES quasar spectra (Q0 109 + 200, 
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TABLE 2 

Galaxy Asymmetries 
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a The D z computations for the galaxies Q0235 + 164 at z = 0.5240 and Q1623 + 268 at z = 0.8875 failed. 



Q0235 + 164, Q0836+ 113, and Q1332 + 552), and cannot di- 
rectly compute their f sat . Fortunately, the Mg II A2796 profiles 
of the absorbers in the Q023 5 + 164 and Q0836 + 113 fields 
are available for inspe ction in lLanzetta & Bowenl (11992b and 
Turnshek et al.l (119891) . respectively. Their spectra are of com- 
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Fig. 2. — The Mgn A2796 equivalent width as a function of the ratio 
of the number of saturated pixels to total pixels, f sat , for profiles observed 
with HIRES and UVES spectra. Both the dashed lines were obtained by 
maximizing the number of points in the lower left and upper right quadrants 
(Ql and Q3) and minimizing the number of points in lower right and upper 
left quadrants (Q2 and Q4). We define a profile as highly saturated for f sa t > 
0.5. The horizontal dashed line provides a natural and objective equivalent 
width transition of W r (2196) = 1.39 A for defining our subsample. 



parable resolution to HIRES and UVES. These two absorbers 
appear to clearly have f sat > 0.5, indicating that they are also 
highly saturated. Based upon our investigation, it is reason- 
able to assume that the very large W r (2196) absorbers in the 
Q0109+200 and Q1332+552 fields also have highly saturated 
profiles. 

4.2. Examining a Homogeneous Subsample 

Applying the above motivated 50% saturation limit, which 
corrersponds to W r (2196) < 1.39 A, we obtain a subsample 
of galaxy-absorber pairs with a higher degree of homogene- 
ity that can be investigated in more detail for a deeper un- 
derstanding of the nature of A/D- W r (2196) correlation. The 
subsample comprises 28 absorber-galaxy pairs from which 
the highly saturated profiles have been excluded. For these 
highly saturated profiles, the single valued equivalent width is 
a less reliable measure of the quantity and complex distribu- 
tion of absorbing gas. 

As with the full sample, we performed non-parametric 
Spearman and Kendall correlation tests with the computed 
asymmetries for the subsample of 28. In addition to the 
equivalent widths, we included high resolution kinematic 
properties of the absorption, i.e., the velocity spread, ve- 
locity asymmetry, number of clouds (based upon Voigt pro- 
file decompositio n) , and the tot al MgH column density [see 
IChurchill & Vo^tl d200lh and IChurchill. Voet. & Charlton] 
(120031) for definitions of these quantities] . Though these kine- 
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matic quantities correlate with 1^(2796) to varying degrees 
(only the number of clouds is statistically significant), no sin- 
gle kinematic quantity provides a unique description of the 
magnitude of absorption as measured by the equivalent width 
dChurchill & V ogt 2001). Interestingly, none of the kinematic 
absorption properties were found to correlate with the vari- 
ous asymmetry measures, nor with the asymmetries normal- 
ized by impact parameter. However, the significance of the 
A/D- W r (2196) correlation remains highly significant with a 
confidence level (CL) of 99.9% (a 3.33 a result that the null- 
hypothesis of no correlation is ruled out). 

We explored the statistical behavior of the correlations be- 
tween WX2796) and all the asymmetry measures computed 
using GIM2D. Recall that the Abraham A is computed over an 
elliptical aperture defined by 2 a above background, whereas 
Ra and A z are computed for circular apertures at 10 values of 
Vh (model dependent) and for three different levels of signif- 
icance. A selected subset of our exploration is presented in 
Table [3] In column (1) are the tested properties. In columns 
(2)-(4) are the Spearman correlation coefficient, r s , the prob- 
ability, Ps, that the tested data are consistent with the null hy- 
pothesis, and the number of standard deviations, N a assuming 
the ranks are drawn from a normal distribution. In columns 
(5)-(7) are the Kendall r, the probability, P%, and N a . For 
simplicity, we will hereafter discuss the Kendall statistics. 

In Table [3l we first report that there is no significant corre- 
lation (or anti-correlation) between impact parameter, D, and 
W r (2796) for our subsample. The equivalent width does trend 
toward smaller values as D increases, but the CL is 92.5% 
(1.8 a), which does not classify as a statistically significant 
departure from the null-hypothesis. However, in our explo- 
ration, we consider the possible affects of the trend between 
D and W r (2796). 

We first focus on the asymmetry Ra for both R 3 A ay2r ' 2 and 
£3a,3r h ^ ^ a b ove background for 2 and 3 respectively). As 
presented in Table there is a total lack of a correlation for 
each of W r (2796) versus R 3 / M and W r (2796) versus R 3 / M 
with r K = 0.07 and 0.06, respectively. These statistics trans- 
late to confidence levels (CL) of 40.2% and 34.5%. When the 
Ra are normalized by D, moderately significant correlations 
arise at the 98.7% (2.5 a) and 99.3% (2.7 a) CL, respectively. 
To investigate these results as a function of the number of i%, 
we present the behavior of the CL for both R A a /D and R A a /D 
as a function of the number of half-light radii out to 5 
in Figure [3]2. For both sensitivity thresholds, the CL scatter 
around ~ 99%, which are moderately significant correlations. 
The CL for the Ra substantially decrease for 1 r^. This is not 
unexpected since the Ra are computed using the model resid- 
ual images; the GIM2D models have difficulties matching the 
surface brightness profiles in the galaxy cores, resulting in ar- 
tifacts in the residual images for r& < 1 . These model depen- 
dent artifacts tend to randomize the Ra in the cores, which 
washes out any signal in the correlation tests. 

We next list a selected correlation test for the asymmetry for 
^2aMh (pixels 2 a greater than their symmetric counterparts 
within \rh) and the correlation test for D z in Table [3] As with 
Ra, when the A z is normalized by D, the Kendall probability 
increases from a statistically insignificant CL of 79.6% (1.3 a) 
to a marginally significant 96.8% (2.1 a). In Figure [3?, we 
show the CL behavior of A 2 z a /D and A 3 z a /D as a function of 
the number of half-light radii out to 5 . The A z are computed 
directly on the data; there are no artifacts at < 1 introduced 
by the GIM2D models. The CL of the A 2 z a /D and A^/D 



TABLE 3 

Results of K-S Tests With W r (2196) < 1.39 A 



Spearman Kendall 
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K 


W r vs. D 


-0.31 0.1137 1.58 -0.24 


0.0746 


1.78 


Ra Method (Schade et al. 1995) 


W r VS. R 3 /^ 


0.06 0.6798 0.41 0.07 


0.5977 


0.53 


W r VS. flf 7 '^ 


0.09 0.5713 0.57 0.06 


0.6552 


0.45 


W r vs. R 3 / arh /D 


0.45 0.0181 2.36 0.33 


0.0126 


2.50 


W r vs. R 3 / M /D 


0.50 0.0095 2.59 0.36 


0.0071 


2.69 




and D 7 Methods (Simard et al. 2002) 






W r VS. A z h 
W r VS. D z 


0.23 0.2082 1.26 0.17 
0.03 0.8319 0.21 0.02 


0.2041 
0.8885 


1.27 
0.14 


W r YS.A 2 z aMh /D 
W r VS. D z /D 


0.43 0.0243 2.25 0.29 
0.13 0.5029 0.67 0.08 


0.0324 
0.5523 


2.14 
0.59 


A Method (Abraham et al. 1994) 


W r vs. A 


0.40 0.0359 2.10 0.28 


0.0394 


2.06 


D vs. A 


0.25 0.1889 1.31 0.20 


0.1437 


1.46 


WrVS.A/D 


0.64 0.0008 3.36 0.45 


0.0009 


3.33 



tests substantially decrease for increasing r^. The correlation 
test for D z and W r (2196) reveals no indication of even a slight 
trend (yielding r K ~ 0.02) even when D z is normalized by D 
(yielding r K ~ 0.08). Recall that D z is sensitive to the presence 
of tidal arms for > 1 . 

The different behavior of the CL for the Ra and A z cor- 
relation tests can be explained by considering the effect of 
Poisson noise. Ra is sensitive to significant fluctuations about 
a smooth, correlated, axisymmetric model surface. As the 
apertures increase, the Poisson fluctuations in the pixel values 
become relatively more important, but at a rate much slower 
than for that of A z . This is because A z is based upon differ- 
ences between pixel values and the Poisson fluctuations at 
larger r h apertures can wash out the significance of these dif- 
ferences. 



-(a) ' 


Rf/D : 

Rf/D - 























Af/D 



Half Light Radius, r h 



Half Light Radius, r h 



Fig. 3. — (a) Confidence level (CL) for the correlation between Wr(2796) 
and Ra/D over the range 1 < < 5 for 1 a (solid line) and 3 a (dashed 
line) thresholds above background. R& is computed from the model residual 
images, which are less robust near the galaxy core (r^ < 1). — (b) CL for 
the correlation between W r (2796) and A z /D over the range 1 < r/, < 5 for the 
2 a (solid line) and 3 a (dashed line) cases. A z is computed directly from the 
data. Both R A and A z employ the circular r/ 7 apertures of the models. Note 
the different CL scales for panels a and b. 

The correlation test for the Abraham A with W r (2796) is 
also listed in Table [3J There is a moderately significant cor- 
relation between W r (2196) and A, with CL of 96% (2.1 a). 
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W r (2796) A 

Fig. 4. — Galaxy morphological asymmetry A normalized to impact pa- 
rameter, D, versus W r (2196) for the subsample with f sat < 0.5. The error 
bars are 1 a; in most cases, the errors are smaller than the data points. The 
correlation has a 99.92% (3.3 a) confidence level. The dashed line is a max- 
imum likelihood fit, A/D = 2.483 x 1 0~ 3 W r (2796) + 8.15 x 1(T 4 . — (inset) 
The same plot for the full sample, illustrating the increased scatter for larger 
equivalent widths. The dashed line is the cutoff for the subsample. 

We find only a weak trend between A and D (1.5 a). How- 
ever, when the A are normalized by D, a significant correla- 
tion emerges with a CL of 99.9% (3.3 a). Since the A are 
computed within a single isophotal area, we cannot examine 
the behavior of A with r^. In Figure |U we present A /D versus 
WX2796), with the dashed line showing a maximum likeli- 
hood fit. 

It could be argued that the correlation between A/D and 
Wr(2796) might be induced by the weak trend (1.8 a) be- 
tween WX2796) and D. To investigate this, we selected corre- 
lation tests of other GIM2D model galaxy parameters versus 
Wr(2796) for which the CL from the correlation tests were 
similar to the tests of the asymmetries presented in Table [3] 
We then divided these model galaxy parameters by D and 
reran the correlation tests. None of the scaling tests yielded an 
"induced" correlation (the highest significance was ~ 2.0 a). 
In fact, in all cases where 2.0 < N a < 2.4 prior to normalizing 
by D, these scaling tests yielded reduced significance. 

To further investigate the confidence level we found for the 
A/D-W r (2796) correlation, we applied the bootstrap replace- 
ment program resampling dHowell 2004) to the subsample. 
We performed 10,000 bootstrapped samples, each of n = 28, 
with random replacement from the original data to determine 
if the correlation is due to a chance realization of the data. The 
bootstrap statistic corroborates the Spearman-Kendall statis- 
tics, yielding a 99% confidence level. 

5. DISCUSSION 

In our exploration of the morphological properties of galax- 
ies selected by Mg II absorption, we have found a statistically 
significant correlation (3.3 a) between the Mgll A2796 equiv- 
alent width and galaxy asymmetry, A, normalized by the im- 
pact parameter. Background fluctuations in the calculation of 
A introduce only a 0.1 a spread in the significance level of 
the correlation. Additional uncertainties in the A values due 
to PSF subtraction of the quasar in the seve n small impact pa- 
rameter galaxy-absorber pairs (see § 13.11) introduced only a 
0.2 a spread in the significance level. 

It is of interest to address the individual correlation test re- 
sults for A versus D and for A versus W r (2196) as presented in 
Table [3] First, A versus D shows a weak positive trend such 
that larger A are found at larger D. This suggests that uncer- 
tainty due to PSF subtraction of the quasar is not inducing 



elevated A values for galaxies at small D. Second, A versus 
Wr(2796) exhibits a moderately significant correlation (2.1 a) 
in that large W r (2796) are associated with galaxies measured 
to have larger A. This hints that there may be a direct connec- 
tion between galaxy morphological asymmetry and the quan- 
tity of gas in the galaxy halo. 

However, as described above, normalization of A by D ac- 
centuates that trend, yielding a significant correlation. There- 
fore, for absorption selected galaxies, the A/D- W r (2196) 
correlation suggests a connection between the magnitude of 
asymmetric perturbations in a galaxy's morphology, normal- 
ized to the galactocentric distance probed by the quasar, and 
the quantity of absorbing gas residing at that distance in its 
halo. In effect, this might be interpreted such that the in- 
fluence of galaxy asymmetry diminishes with distance in a 
halo. The significance of the correlation is substantiated by 
the bootstrap replacement method. We note that this result is 
independent of the GIM2D galaxy models, since A is com- 
puted directly from the counts in the galaxy. 

5.1. Comparison of Asymmetry Measures 

The different computational details of A, Ra, A z , and D z 
render each with varying degrees of sensitivity to the details 
of the underlying structures that give rise to the morpholog- 
ical asymmetries in the galaxies. Furthermore, the relative 
behavior between the correlation tests provide insight for in- 
terpreting the correlation. 

The CL behavior of A z with half-light radius aperture, r&, 
suggests that the tendency for a correlation of A z /D with 
Wr(2796) arises in the centrally concentrated region of the 
galaxies, primarily within 1 (see Figure 0. In contrast, 
the CL behavior of Ra is fairly flat as a function of the r& cir- 
cular aperture. For larger apertures, the sensitivity of Ra to 
asymmetries (based upon model residuals and the parameter 
ru) may be higher than that of A z (based upon the data and the 
parameter r&). However, since Ra is an entirely model depen- 
dent quantity, it is not highly reliable in the central aperture 
(jh = 1) due to difficulties with the model in the galaxy cores 
(yielding poor subtraction of the galaxy model in the core and 
thus resulting in a low CL). The low CL of D z , which is com- 
puted exclusively outside the r^-X, suggests that the asym- 
metries are not due to strong spiral arms or tidal arms to one 
side of the galaxy (which is both consistent with the behavior 
of A z and visual inspection of the galaxy images). 

On the other hand, the Abraham A is an entirely model in- 
dependent quantity. It is computed directly on the data over 
a generalized elliptical aperture defined by the region where 
galaxy counts are 2 a above background. Thus, the A/D- 
Wr(2796) correlation is independent of any arbitrary assump- 
tion (such as the functional form of the bulge and disk surface 
brightness profiles) applied in the GIM2D modeling. Since 
the apertures are operationally defined using the signal-to- 
noise ratio of the galaxies, the asymmetries could in princi- 
ple be sensitive to the pixel counts of each galaxy. The pixel 
counts for a given galaxy depend upon the product of expo- 
sure time and surface brightness. We calculated the signal- 
to-noise ratio from the surface brightness of each galaxy over 
their full 1.5 a aperture. A correlation test yielded no trend 
between A and galaxy signal-to-noise ratio. 

The combined behavior of the asymmetries and the corre- 
lation tests with rh indicate that the asymmetries of the galax- 
ies in our sample are centrally concentrated. The asymme- 
tries are not the result of major perturbations that would result 
in significant departures from a de Vaucouleurs profile in the 
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bulge components and an exponential profile in the disk com- 
ponents. As such, the morphologies of Mgll absorption se- 
lected galaxies appear to be "normal" in their classifications. 
This suggests that extraordinary events, such as major galaxy 
merging or a history of galaxy harassment, are not character- 
istic of galaxies hosting extended gaseous halos. 

5.2. Asymmetries and Absorption Selection Galaxies 

Since the galaxies in our study are selected by association 
with Mg II absorption, it is of interest to know if this selection 
criteria results in a population of galaxies with higher levels of 
morphological asymmetries (or perturbations) than the gen- 
eral population of magnitude selected galaxies at intermediate 
redshifts. The Medium Deep Survey (MDS, lAbraham et al.l 
1996b), which employs the F814W filter, is the only survey 
to which we can undertake a fair comparison since the MDS 
galaxies have also been modeled using GIM2D and are drawn 
from both a similar magnitude range (20 < mFxuw < 2 2) and 
redshift range as our sample. ICrampton et al.l (Il995|) argue 
that the majority of galaxies in the MDS reside at z < 1.0 with 
a peak in the distribution at z — 0.6. The range of magnitudes 
of our sample is 19.4 < mp^w < 23.8 and covers the redshift 
interval 0.3 < z < 1 with a peak in the distribution at z = 0.43. 

It is well known that the details of galaxy morphologies 
depend upon the rest-frame bandwidth in which they are ob- 
served. Almost all of the galaxies in our sample were ob- 
served with the F702W filter; only two have been observed 
in the F814W filter, whereas the MDS galaxies are all ob- 
served in the F814W filter. Thus, there is some concern that 
a quantitative comparison of asymmetries between our sam- 
ple galaxies with those of the MDS may be compromised 
due to a non-uniform rest-frame band in which the galaxies 
are observed. We can only address this concern statistically. 
Since the peak of the redshift distribution of MDS galaxies 
is z ~ 0.6, the majority of the galaxies are observed in the 
rest-frame wavelength band 4530-5490 A. For our sample, 
the redshift distribution peaks at z = 0.43, which translates to 
a rest-frame wavelength band of 4360-5340 A for the major- 
ity of our sample. Therefore, there is a significant rest-frame 
wavelength overlap between the two samples. 

Given these caveats, we performed a K-S test to examine 
whether the Abraham A asymmetries for the galaxies in our 
sample are consistent with having been drawn from the distri- 
bution resulting from the MDS survey. In Figure [3 we show 
the normalized distribution of A. The mode of the asymme- 
try distribution for our sample in = 37) is highly peaked at 
A = 0.08 whereas the mode for the MDS sample (n = 435) is 
A = 0.04. The K-S probability that the two samples are drawn 
from the same A distributions is P KS = 0.0023, which rules out 
the null hypothesis of similar distributions at the 99.77% con- 
fidence level (2.8 a). 

This is a remarkable result given the uncertainties in the di- 
rect comparison. The dominant uncertainty is due to our lack 
of kn owledge of the redshift s of the galaxies in the MDS sam- 
ple. [Abraham et alJ (fl996b) suggest that the extended tail in 
their A distribution is dominated by z > 1 galaxies, which are 
observed in the rest-frame U-band where morphological per- 
turbations are significantly pronounced. Without redshifts for 
the MDS galaxies, we cannot select out the z > 1 galaxies for 
a more robust test. However, it is clear that the high redshift, 
high asymmetry tail of the MDS galaxy distribution results in 
a lower confidence level for the K-S test than if the tail could 
be objectively corrected for redshifting effects. 




Fig. 5. — The normalized distribution of asymmetries (Abraham's A) for 
the MDS (dotted line) and our sample (solid line). 

It is important to note that there only two Mg II absorption 
selected galaxies with A < 0.04, below the mode of the MDS 
distribution. This constitutes only 5% of our sample. How- 
ever, 29% of the A distribution for the MDS sample is popu- 
lated below this value. It is the difference in this portion of the 
cumulative distribution function of A that dominates the K-S 
statistics and strongly suggests that Mg II absorption selected 
galaxies tend to have larger morphological asymmetries than 
those of the galaxies in the MDS. It is quite possible that there 
is a bias towards galaxies with slightly larger asymmetries us- 
ing Mgll absorption selection. 

5.3. Interpretation and Context of the Correlation 

Normalizing a galaxy property by the projected distance 
to the quasar provides a measure of the influence the galaxy 
property has out to the extended regions where the gas is by 
chance sampled. If the influence falls off with distance, then 
normalizing by D provides a "scale length correction". Our 
results imply that galaxy asymmetry is an important quan- 
tity for predicting the abundance of halo gas and its covering 
fraction and that the influence of asymmetry strengthens with 
proximity to the galaxy. 

For local galaxies there is evidence that morphologi- 
cal perturbations are associated with a high covering frac- 
tion of gas beyond the v i sible components of the galax- 
ies (e .g. JPuche et all 1 19921: ISwaters. Sancisi & van der Hulst 



1997; Rand 2000; Fraternali et al. 2001, 2002; Barbieri et al 



2005, also see Braun 2005). Our results for intermediate red- 
shifts are consistent with this trend. Apparently, the greater 
the morphological perturbations of a galaxy, the more efficient 
the process of distributing gas into its halo. Furthermore, our 
results suggest that similar absorption strengths can arise at 
larger galactocentric distances from galaxies with more pro- 
nounced perturbations or at smaller galactocentric distances 
from galaxies with less pronounced perturbations. 

Perturbations at the level measured here can be produced 
via minor satellite galaxy mergers with the host absorbing 
galaxy, mild (long range) galaxy harassment, or continued dy- 
namical interaction with (i.e., accretion of) the local cosmic 
web. Satellite galaxy interactions can generate Magellanic- 
type streams, increasing the gas covering factor out to ~ 50 
kpc. These various mechanisms likely also gently modify the 
galaxy morpholog y and induce stochast ic periods of elevated 
star formation (Kauff mann et al.l l2006|) . These episodes of 
star formation generate galacti c fountains, ch imneys, super- 
bubbles, and low-level winds (Dahlem 1997). Galaxies ex- 
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hibiting negligible morphological perturbations may not have 
experienced these processes within the timescale that gas cy- 
cles in halos. These galaxies would be expected to produce 
very weak or undetectable Mg II absorption even at moderate 
impact parameters since they would be deficient in mecha- 
nisms that populate their halos with gas. 

This scenario would explain why some bright galaxies clos e 
to the quasar produce weak absorption ( Churc hill et al.ll2 006). 
A natural consequence of this scenario is that it also predicts 
some galaxies might have such low quantities of gas in their 
halos that they would be excluded from an absorption selected 
sample even if they lie at relatively small impact parameter to 
a quasar. In the course of our studies, we discovered a galaxy 
at z = 0.6720 in the Q1317 + 274 (TON153) field at an impact 
parameter of D = 58 kpc for which the 3 a upper limit on ab- 
sorption is W r (2796) = 0.007 A (iChurchill et alJl2007h . The 
galaxy asymmetry is A = 0.08, giving A/D = 0.0014. This 
data point, though an upper limit, is consistent with expec- 
tations from the A/D- W r (2196) correlation, given the scatter 
of the data about the maximum likelihood fit. If the A/D- 
WX2796) correlation holds true, it would be expected that 
"non-absorbing" galaxies have very small A/D. 

The galaxy-absorber pairs for which the absorption profiles 
are more than 50% saturated exhibit a higher level of scatter 
in the A/D- W r (2196) correlation than do those that are not 
highly saturated. This may in part be due to the loss of infor- 
mation in the combined velocity spread and column density 
distribution of the gas in the highly saturated profiles. Our 
objective examination of profile saturation resulted in a natu- 
ral equivalent width break of WX2796) ~ 1.4 A. This equiva- 
lent width regime may mark a true transition in the underlying 
physical mechanisms giving rise to the gas. Recently, we have 
found that galaxies selected by "weak" Mgii absorption are 
statistically consistent with those selected by so-called "clas- 
sic" Mgii absorbers (IChurchill. Kacprzak. & Steidell [20051: 
IChurchill etaDl2006h . For larger W r (2796), there is a grow- 
ing body of evidence that DLAs, wind-driven galaxies, and 
strongly interacting galaxies are selected. 

Using a statistical argument, IZibetti et all (l2006|) demon- 
strated that galaxies selected by Mgii absorption with 
Wr(2796) > 1 . 1 A have colors consistent with star burst galax- 
ies, whereas galaxies selected by W r (2796) < 1.1 A are con- 
sistent with redder, less active galaxies. The mean color is 
consistent with that of an Sb spiral. Also using statistical 
methods, iBouche et al.l (l2006|) reported an anti-correlation 
between the mass of Mgii absorbing halos and W r (2196). 
They argue that this correlation is inconsistent with virial- 
ized gas and suggest that a large fraction of strong Mg II ab- 
sorbers galaxies [those with W r (2196) > 2 A] are produce by 
galactic winds and/or outfl ows produc ed by active s t ar for- 
mation. iBond et al.l (1200 ll) argue and iProchter et all (12006b 
similarly suggest that Mgii absorbers with W r (2196) > 1.8 A 
arise in galaxies with elevated star forma tion, and that the gas 
is likely wind driven. iNestor et all (|2006|) imaged the galaxies 
in the fields of 15 quasars with "ultra" strong Mgii absorp- 
tion [2.7 < W r (2796) < 6.0 A] and found that, compared to a 
control sample, there is a significant excess of L > few x L* 
galaxies and also possibly interacting galaxies in the quasar 
fields. 

DLAs are also a very different environment than the 
Lyman-limit and sub-Lyman limit absorbers in our subsam- 
ple. Highly saturated Mgii profiles with modera te velocity 
spreads preferentially select DLAs (Ellison 2006). It is also 



observed ( Rao 2005) that the fraction of Mg II absorbers that 
select DLA environments increases with increasing W r (2196). 
However, it is common for DLAs to preferentially select 
smaller impact parameters, which would systematically ele- 
vate their A/D values. 

Therefore, we can infer that absorbers with f sat < 0.5 arise 
predominately from galaxies that have experienced either 
mild levels of interactions and/or minor mergers or ongoing 
accretion of the cosmic web. These processes would increase 
the level of star formation. The resulting asymmetry may re- 
flect both the gravitational disturbance of the galaxy as well 
as surface brightness variations due to induced star formation. 
Absorbers with f sat > 0.5 may arise predominantly from ei- 
ther DLAs, strong interacting galaxies, or galaxies with strong 
winds. If so, it would not be unexpected that active star form- 
ing galaxies hosting winds might also adhere to the A/D- 
WX2796) correlation, but with a larger scatter due to their 
highly saturated absorption profiles. 

6. CONCLUDING REMARKS 

We have examined four quantified measures of morpholog- 
ical asymmetries, A, Ra, A z , and D z in a sample of 37 Mgii 
absorption selected galaxies. The Mgii absorbers were se- 
lected from a database of HIRES and UVES quasar spectra, 
and from the literature in a few cases. We modeled the galax- 
ies, observed with WFPC-2/HST in the F702W and F814W 
filters, using GIM2D. The Abraham A is model independent, 
whereas Ra is based upon the residuals of the GIM2D mod- 
els. The A z and D z are computed directly from the data, but 
invoke apertures determined from the GIM2D models. 

We report two main results: 

1. We find a statistically significant correlation between 
the Abraham A normalized by the quasar galaxy impact 
parameter, D, and the Mgii A2796 rest-frame equiva- 
lent width. The significance is 3.1 a for the our full 
sample of 37 galaxy-absorber pairs. However, when 
we limit the sample to those in which the fraction of 
saturated pixes in the Mgii A2796 absorption profiles 
is fsat < 0.5, the A/D-W r (2196) correlation is signifi- 
cant at the 3.3 a level. 

2. A Kolomorov-Smirnov test comparing the distribu- 
tions of Abraham A for Mgii absorption selected galax- 
ies and magnitude selected galaxies from the Medium 
Deep Survey shows the distributions are not consistent 
at a confidence level of 99.8%. There are caveats re- 
lated to the comparison; however, the data show a clear 
paucity of small A for Mg II absorption selected galax- 
ies. Selecting galaxies by absorption in quasar spectra 
may bias against galaxies with negligible morphologi- 
cal perturbations. 

The A/D-W r (2796) correlation is formally statistically sig- 
nificant, if only at the 3.3 a level. However, it is a quite re- 
markable finding given the large vari ations in kinematic prop - 
erties of Mgii absorption profiles (|Churchill & Vo gt 2001). 
As discussed in § [3 the behavior of the correlation is consis- 
tent with a growing body of evidence that the quantity of gas 
in galaxy halos is connected to the processes occurring in the 
galaxies themselves, and that stronger, highly saturated Mgii 
absorbers are related to higher levels of activity. 

A larger sample of galaxy-absorber pairs with both HST 
imaging data and HIRES and/or UVES quasar spectra would 
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be instrumental for confirming and possibly strengthening the 
statistical significance of our results. One of the major pitfalls 
of studying absorption selected galaxies is that it requires in- 
tensive work to obtain a thorough census of the galaxy red- 
shifts in quasar fields. Thus, we cannot claim that the galax- 
ies used in this study are the only galaxies associated with the 
Mgll absorbers. 

Our results suggest a relationship between the gaseous halo 
and the morphological perturbations in the host galaxy nor- 
malized to the impact parameter of the observed absorp- 
tion. Further insight may be obtained by comparing the 
relative kinematics of the absorbing gas with those of the 
host galaxy. Only a half-dozen absorptio n selected galaxies 
have been studied at this level of detail (Ste idel et al.ll2002t 
Ellis on et"aT1 l2003). In addition, further exploration of corre- 
lations between C IV ab sorption strength and Mg II kinematics 
( Churchill et al. 1999a) with the galaxy morphological pertur- 
bations would provide important clues for better understand- 
ing the A/D-W r (2196) correlation. We are currently working 
toward these goals. 
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